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Abstract 

For low T new strict results for the instanton density n(T) are reported. Using 
the PCAC methods, we express n(T) in terms of vacuum average values of certain 
operators, times their calculated T-dependence. At high T, we discuss the applicabil- 
ity limits of the perturbative results. We further speculate about possible behaviour 
of n(T) at T ~ T c . 
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1 Introduction 



Tunneling between topologically different configurations of the gauge field, described 
semiclassically by instantons |l|], dominate the physics of light quarks. In early works 
(summarized e.g. in [0]) instantons were treated as a dilute gas, while later it was 
recognised that the instanton ensemble resemble rather a strongly interacting "in- 
stanton liquid" |5[ [| . During the last years calculations of the correlation functions 
|?J and Bethe-Salpeter wave functions |8| for various mesonic and baryonic channels 
were made along these lines. The results agree surprisingly well both with phe- 
nomenology H and lattice simulations |[Tof1 . Parameters of the "instanton liquid" 
were also reproduced (by the "cooling" method) directly from the (quenched) lat- 
tice configurations (O]. In addition to that, it was found that correlation functions 
as well as hadronic wave functions in most channels remain practically unchanged 
after "cooling". In particular, main mesons and e.g. the nucleon remains bound, 
with about the same mass and wave function. This confirms that the agreement 
of previous lattice calculations with the instanton model was not accidental, and 
instantons indeed are the most important non-perturbative phenomena in QCD. 

Investigations of the finite temperature case were started in [12], where classical 
caloron solution was found. Although the solution depends on T, the action is T- 



independent. Furthermore, it was argued by one of us [14] that at high T the 
specific charge renormalization and the Debye-type screening of the electric field in 
quark-gluon plasma should suppress instantons with size p > 1/T. Pisarski and 
Yaffe |22][ have evaluated the T-dependence in the one-loop approximation. The 
physical nature of their result and its applicability region will be discussed below. 

Last years the studies of the finite temperature had focused especially on the 
region around the chiral phase transition T ~ T c . The first attempt to understand 
this phase transition as a rearrangement of the instanton liquid, going from a ran- 
dom phase at low temperatures to a strongly correlated "molecular" phase at high 
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temperatures was made in [jD?. Recently this idea was recently made more quanti- 
tative in JT||, [l6|], and although the detailed comparison to lattice thermodynamics 
and correlation functions is yet to be made, the first results show overall agreement, 
indicating that the mechanism of chiral restoration is basically understood. 

The particular topic of this paper, the temperature dependence of the instanton 
density n(T), is certainly an important ingredient of all this development. However, 
it has attracted surprisingly little attention in literature. The only attempts to 
determine this quantity by "cooling" of the lattice configurations (same method as 
was used at T=0) were made in refs. JlTj . The only statement one can probably make 
using these data is that the density has no significant T-dependence, till T ~ T c . 
Unfortunately, accuracy of this statement remains at the level of 50%, at best. 

2 Instantons at low temperatures 

As far as we know, the modification of the instanton density in this limit was never 
considered before. However, the general physical picture at low T « T c is well 
known: the heat bath is just a dilute gas of the lightest hadrons, the pions. The 
problem is especially clear in the chiral limit, in which quark and pion masses are 
neglected, and this case is assumed in what follows. Due to their Goldstone nature, 
the large wavelength pions are nothing else but small collective distortions of the 
quark condensate. Therefore, one can always translate the average values over the 
pion state as the vacuum expectation value (VEV) of a different but related quantity. 

Let us now consider a tunneling event, described semiclassically by the instanton 
solution. As discovered by 't Hooft ||, it can only take place if certain rearrange- 
ments in the fermionic sector are made, which can be described by some effective 
Lagrangian with 2Nf fermionic legs. For simplicity, in this work we assume that 
the number of light flavors Nf = 2, disregarding strange and heavier quarks. 

The situation can be further simplified by consideration of small-size instantons 
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P^-QCD « 1) for which this Lagrangian can be considered as a local operator. What 
follows from 't Hooft Lagrangian, after averaging over the instanton orientations is 
made, is [18|: 

AC = J dpd (p)(^ 2 p 3 f{qT + qqT^q + ^qT a + qqT<Lq - ±-qT a + ^ qqT a _^ q} (1) 
where the definition of the operators involved is as follows 

r ± = (^)®(^p) (2) 



n = { l -^)®i l -^)®t a (3) 



rr = (^)<v®(^)®t° w 

At T=0 the instanton density is therefore proportional to the VEV of AC, and 
the only thing which changes at low T is clearly the modification of the quantities 
above. 

How to do this technically was actually clarified by PCAC-related paper in 60's. 
The necessary formulae can be found e.g. in the recent paper by Eletsky |l8|| , where 
different set of four-fermion operators (appearing in QCD sum rules for vector and 
axial currents) was considered. The general expression is 

((qAq)(qBq)) T = {(qAq)(qBq)) - ^{{^1 {r a B A}}q)(qBq)) 



~ ^L((q{nA}q)(q{r$B}q))o (5) 



96F 2\v^ JL J JJ ^ ;/U 48F2 

where A,B are arbitrary flavor-spin-color matrices and T| = T a ^ B . Here and below 
flavor matrices are shown as r a , and color one as t a . 

For generality, there are six different operators of different spin-flavor-color 
structure involved (see the Table). Their T-dependence in 0(T 2 ) order can be found 



from the expression above, and it is also listed in the Table. Generally, the operators 
mix, and it is convenient to group those combinations which do not. Returning to 
the effective Lagrangian, one can see that there are only two combinations which 
are actually relevant 

The final result for the instanton density at low temperature T therefore contains 
two constants, the vacuum averages of these operators 



P 



IT 2 IT 2 



(8) 



Although two VEV's which appear here are unknown (and are subject to further 
investigations), it is clear that the total T-dependence should be rather weak: it is 
bound to be 1 + a-£| with a in the strip a = (—1/6, 1/6). 

The so called vacuum dominance (VD) hypothesis [JTs[] was used in various 
applications (such as QCD sum rules and weak decays) for evaluation of VEV's of 
various operators. It leads to VEV's and the 0(T 2 ) corrections also indicated in 
the Table 1. Remarkably enough in this case the T-dependence exactly cancel. 

Returning to discussion of our general result (8), we comment that this result by 

itself rules out some possible picture of low-T vacuum structure. In particular, the 

so called "random instanton liquid model" (RILM) was shown to be a reasonable 

approximation for the T=0 case 0. One may wander if the same model can describe 

the T-dependence at low T. If it is the case, the quark condensate should scale with 

the instanton density as n]^ t , see ||]. Being combined with the well known chiral 
1 (However, the reader should be warned that in general this approximation is not supposed 
to hold or be very accurate, particularly for the operators considered, which are related with 
instantons. 
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theory result 

(qq)T = (gg)o(l - -^) (9) 

these two formulae lead toQ rii ns t (T)/n inst (0) = (1 - 

However, as this estimate happen to be outside of the strip indicated above, this 
possibility is definitely ruled out. It means that, even if RILM is a perfect model 
at T=0, it cannot be so for even small T. This conclusion agrees very well with 



other studies of the instanton ensemble, such as [15|, which emphasize the role of 
correlations built up with growing T in the ensemble of instantons. 



3 Instantons at high temperatures 



QCD vacuum at high temperatures undergoes a phase transition into a new phase, 
called the quark-gluon plasmaffij^ . Although virtual gluons antiscreen the external 
charge (the asymptotic freedom), the real gluons of the perturbative heat bath 



screen it, leading to the well known expression for the Debye screening mass [ 20 1 : 



M 2 D = (N c /3 + N f /6)g 2 T 2 



(10) 



where N c ,Nf are the numbers of colors and flavors, respectively. "Normal" 0(1) 
electric fields are therefore screened at distances 1/gT, while stronger "non-perturbative" 
fields of the instantons 0(1/ g) should be screened already at scale 1/T. 

Quantitative behaviour of the instanton density at high temperatures was de- 



termined in ref. [22 



dn(p, T) = dn(p, T = 0) exp \ -\\ 2 (2N C + JV» - 12.4(A) 



l + l(N c -N f ) 



(11) 



2 Note that the same result can be obtained by a naive assumption, which was used in some 
works on QCD sum rules in the past: namely, that average of all four-fermion operators have 
T-dependence as the square of the condensates. 



6 



where A = vrpT, and A(X) = — ^A 2 + o(A 2 ). Therefore, at high temperatures the 
contribution of small size instantons such as T >> 1/p is exponentially suppressed. 
As a result, the instanton- induced contribution to physical quantities like energy 
density (or pressure, etc) become of the order o|] 

e m ~ [ 1/T ^(pA)(H^/3) „ T 4 (A/T) (117V c /3) (12) 
JO P 5 

which is small compared to that of ideal gas e(T)j^ ea ; ~ T 4 . 

Although the Pisarski-Yaffe formula contains only the dimensionless parameter 
A, its applicability is limited by two separate conditions: 

p«l/A, T » A (13) 

The former condition ensure semiclassical treatment of the tunneling, while the 
latter is needed to justify perturbative treatment of the heat bath. In this section we 
would like to discuss applicability conditions of these well-known results in greater 
details. 

Our first point is that the one-loop effective action discussed by Pisarski and 
Yaffe actually consists of two parts with very different physical origin and interpre- 
tation. To show that in the simplest case, consider the determinant corresponding 
to scalar isospin 1/2 fieldQ and rewrite them as follows: 

4 _ r.N ^fl )] - Tr^l-^)] = Sl + S2 (14) 



^^wzW^^^zW), (15) 



3 Here we consider only the pure glue theory. In the theory with massless fermions individ- 
ual instantons are impossible, and only "instanton-antiinstanton" molecules can appear at high 
temperatures. 

4 The determinants of the actual quadratic fluctuations of the quark and gluon fields (modulo 
the factor corresponding to zero fermion modes) can be expressed via the determinants of scalar 
fields with isospins 1/2 and 1 
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^^W^Wd-^W^W), (16) 

Here Tr? is a trace over all matrix structures, plus integration over A4 - the strip 
in R 4 with span in the r direction of 1/T. A(p, T) is the caloron field and A(p) is 
the instanton field. Two contributions introduced in this way, <5i,<52, are the origin 
of two terms in the resulting formula (11). 



As it was shown in ref . 21 ] , the first term can be expressed via the forward scat- 
tering amplitude of heath bath constituents, on the instanton field. Therefore its 
physical origin is clear: p 2 comes from the scattering amplitude, while the temper- 
ature factor T 2 enters via the standard thermal integral over the particle momenta: 

61 = I WT 3 MeMp/T)-l) TrT{P ' P) 
Let us show how it works using the example of a "scalar quark", which is 



simpler than realistic spinor and vector particles considered in [21]. One can evaluate 
T(p,p), the forward scattering amplitude of a scalar quark on the instanton field, 
using standard Leman-Simansik-Zimmermann reduction formula: 

TrT(p,p) = J d A xd 4 y e ip ^ x - y) Tr(d 2 x Ai(x,y)d 2 ) (18) 

where Ai is the (isospin 1/2) scalar quark propagator (23]: 
2 — 

Ai (x, y) = xV+A.ry.rt ^ (ig) 

~ 4tt 2 (x — y) 2 x 2 y 2 (l + p 2 /x 2 ) 2(1 + p 2 /y 2 ) 2 

By rescaling (18) as £ = px, rj = py, subtracting the trace of the free propagator 

and going to the physical pole p 2 = 0, one gets: 

TrT(p,p) = [ ttf^M—?— (Z*_l_l\ = _ 47r 2 p2 (2Q) 



27r 2 (£-r ? ) 2 VfV 2£ 2 2rj 2 , 

As it is just constant, there is no problem with its analytic continuation to small 

Minkowski momenta of scattered quarks, and plugging (20) into (17) we have: 

1 f 1 for periodic fields ~) 

£i = ^A 2 , 77= (21) 
{ —1/2 for antiperiodic fields J 



Note also, that this scattering amplitude has the same origin (and the same depen- 
dence on N c , Nf) as the Debye mass. 

Although formally any result obtained by the perturbative expansion demand 
smallness of the effective charge g(T) « 1, it is not clear in practice what this con- 
dition actually imply. However, we conjecture that accuracy of calculation sketched 
above is controlled by the same effects as the accuracy of perturbative calculation of 
the Debye mass by itself. If so, one can use available lattice studies of the screen- 
ing phenomena (e.g.[p5|1) and check at which T their results start to agree with 
the perturbative formula (10). We then conclude, based on available lattice data, 
that Debye mass and instanton suppression formula (11) should be valid above 
T > T pert = 3T C rs 500 MeV). 

How strong can this suppression be, at that point? Using a canonical "instanton 
liquid" size of the instanton p ~ 1/3 fm, one gets suppression on the level 10 -3 , 
from the 5\ term alone. It suggests a very dramatic behaviour in the interval from 
1 to 3 T c . 

Let us further speculate about the magnitude of 5i contribution for lower tem- 
peratures. In the interval between T c and T per t it is expected on general ground (and 
observe on the lattice) that the Debye mass Mq — > at the critical temperature: 
screening is gone together with the plasma. However, another suppression mecha- 
nism should substitute it below T c , namely the one due to scattering of hadrons on 
the instanton. This is what we have done above for the low-T case, in which only 
the soft pions should be included. At this time, we do not know how to estimate 
this effect including other hadrons. 



Now we turn to discussion of the second term 62 in (14), A(X) in 22], which 
was actually first obtained by Brown and Creamer in p3| . At small A it leads to 
the following correction 

,5 2 = -(l/36)A 2 + (A 2 ) (22) 
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and thus it has the same sign as 5± and parametric magnitude, just numerically 
smaller coefficient. (In the isospin 1 case <5i = (4/3)A 2 and 82 = — (4/9)A 2 + o(A 2 ).) 

The splitting of the variation of the effective action into two physically different 
contributions is a generic phenomenon. This second term has different physical 
origin, because it is connected to a quantum correction to the colored current, times 
the T-dependent variation of the instanton field, the difference between the caloron 
and the instanton. 

Thus, the finite T effects not only lead to appearance of a usual (perturbative) 
heat bath, but they also modify strong {0{\/g)) classical gauge field of the instanton. 
In Matsubara formalism this is described by the non-linear "interference" of the 
instanton field with its "mirror images" , in the (imaginary) time direction. 

Let us conjecture, that for T < T c this suppression mechanism is actually 
irrelevant, for the following reason. It is well known that in this T domain all 
gluonic correlators decay strongly with distance, because all physical "glueballs" 
states are very heavy. It should make any interaction with the "mirror images" (at 
distance (3 = 1/T) virtually impossible. Estimating this effect as exp(—M g [ ue b a u/T), 
where M g i ue i, a u ~ l.QGeV is the mass of the lightest glueball, one gets even at 
T = T c ~ 140MeV a suppression factor of the order of 10 -5 . 

4 Summary and discussion 

We have studied the change of the instanton density at low and high temperatures. 

In the former case, T << T c , we have considered the heat bath as being made 
of dilute soft pions. Applying PCAC methods (in the chiral limit) we have derived 
strict result (8) for the instanton density n(T) at low T. It implies very weak T 
dependence, which agrees with available lattice measurements inside their (so far 
rather poor) accuracy. It also contradicts to some naive models, for example it 
shows that the "random instanton liquid model", presumably a good description of 
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the QCD vacuum, can not be true even at low T. Fortunately, it perfectly agrees 



with the current ideas about the finite-T QCD [15], pointing out that quark-induced 
instanton- antiinstanton correlations are building up with T, till only instanton - 
antiinstanton "molecules" remain for T > T c . 

Our discussion of the high temperatures can be summarized as follows. For very 



high T > T per t ~ 3T C the perturbative result of Pisarski and Yaffe [22] holds. 

Furthermore, we have pointed out that it consists of two parts, 61,62, with dif- 
ferent underlying physics. The first one is directly connected to occupation densities 
of quarks and gluons from the plasma. It is the same effect as lead e.g. to the Debye 
screening mass, so one knows from lattice data at which T this part of the instanton 
suppression can be trusted. We also claimed that it weakens toward T c , but at the 
same time scattering of hadrons on instantons should appear at T < T c , and we do 
not know how to take it into account (except for soft pions). 

The second term 62 originates from the T- dependent variation of the classical 
field, coupled to a quantum correction to the colored current. We expect this effect 
to become exponentially small for T < T c , but we do not know its T-dependence in 
the strip 1-3 T c , where it can be very strong. 

Finally, let us repeat once more, that understanding of the temperature depen- 
dence of the instanton density is of crucial importance for understanding of non- 
perturbative phenomena at and around T c . Surprisingly little efforts has been made 
to clarify this question. In particularly, we call upon lattice community to make 
quantitative measurements of n(T), which can be done by well known methods. 
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5 Table 



Operator 


Coeff. 
in C 


T-rcnormalization 


Vacuum Dominance 
Hypothesis 


Of = 


1 

4 


(Of)o- 


<™>°(132 360^) 


qk^qqk^q 








o£ = 


1 

12 




{ r£{ 36 360^) 


ql^ T a qq :±-To T a q 








o* = 


3 
128 




^#(-64 + 384^) 


q^fqq^fq 




(3{O?) + (O?) )£* 




o* = 


1 
128 


(O^)o- 


« (192 + 384^) 


qlz25 t i T aqqk^L t i T a q 








0? = 


9 
512 




W(768 4608^) 






(3^)0 + ^)0)^ 




0° = 


9 

1536 


(O 2 7 )o- 


<™>°(2304 4608^) 


q X -^a^r a qq^a^r a q 




(3<0?> o + <(9 2 c >o)-g. 
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